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We investigate the dynamics of pressure driven transient flows of incompressible New- 
tonian fluids through circular microtubes having thin elastic walls under the long- 
wavelength and small deformation assumptions, which are valid for many industrial 
and biological processes. An analytical solution of the coupled fluid and solid equations 
is found using Navier slip boundary conditions and is shown to include some existing 
Womersley solutions as limiting cases. The effect of the slip length at the fluid-solid 
interface is analyzed for oscillatory pressure gradients using a range of slip ratio and fre- 
quency parameters. The solutions for elastic and rigid walls are compared for the cases 
with and without slip boundary conditions for a broad range of the relevant parameters. 
It is shown that the elastic behavior of the microtube couples nonlinearly with the slip 
velocity, which greatly enhances the achievable flow rate and pumping efficacy compared 
to the inelastic case. In addition, it is observed that increasing the slip length produces 
less shear stress, which is consistent with the nearly frictionless interfaces observed in 
many microscale experiments. 
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1. Introduction 

Understanding internal fluid flows at small scales is currently a subject of great 
interest for the development of biological and engineering devices and systems at 
the micro- and nanoscales [Beebe et al., 2002; Squires and Quake, 2005]. Since the 
flow dynamics determine the characteristics of the systems, the study of micro and 
nanoflows is of vital importance to developing better technologies. The behavior of 
pulsatile flows at the micro- and nanoscales, in particular, is a largely unexplored 
area [Hansen and Ottesen, 2006; Vedel et al, 2010], and a better understanding of 
such flows is needed for systems with these characteristics. 

Studies of fluid behavior at the micro- and nanoscales has resulted in the growing 
understanding that fluid mechanics at such small scales differs in fundamental ways 
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from fluid mechanics in the macroscopic world [Gad-el Hak, 1999]. One of the fun- 
damental assumptions in fluid mechanics is the no-slip boundary condition: that the 
tangential component of the fluid velocity equals that of the solid at the wall. Re- 
cent advances in micro- and nanoscale experiments and molecular simulations have 
shown that slip of the fluid on the solid surface occurs at small scales, and thus the 
traditional no-slip boundary condition cannot be applied at these scales [Barrat and 
Bocquet, 1999; Baudry et ai, 2001; Choi et ai, 2003; Joseph and Tabchng, 2005; 
Thompson and Robbins, 1990; Thompson and Troian, 1997]. Many other experi- 
mental and theoretical studies demonstrating the slip phenomenon have been re- 
cently reviewed by Whitby and Quirke [2007] and Rothstein [2010], along with other 
small scale fluid dynamics phenomena such as surface roughness, non-wettability, 
dissolved molecules, corrugation and altered shear rate. Neto et al. [2005] explained 
the most evident mechanisms by which both true slip (where fluid molecules are 
effectively sliding on the solid surface) and apparent slip (where the slip occurs 
not at the solid/fluid interface but in the thin depleted fluid/fluid interface layer 
near the solid boundary). For water over hydrophobic surfaces, or for any liquid for 
which the surface is less attractive than the bulk of the fluid (i.e., a non- wetting 
fluid-solid interface) , there is an especially large slip at the fluid-solid interface with 
corresponding drag reduction properties, that is promising for a wide range of fluid 
dynamics and heat transfer applications such as liquid cooling for microelectronics 
[Ou et al., 2004; Rothstein 2010, Watanabe et al. 1999]. 

Whether a fluid is a liquid or gas the magnitude of the slip at the wall can be 
quantifled by the slip length, Ig, at solid surfaces [Thompson and Troian, 1997]. 
Nearly two centuries ago Navier proposed a generalized boundary condition that 
incorporates the possibility of fluid slip at the solid surface. It assumes that the slip 
velocity is proportional to the shear stress at the wall via the slip length. Its equiv- 
alent form for gases agrees with the kinetic theory of gases, which equates the slip 
length with the mean free path of the gas. One of the most important dimension- 
less physical parameters for fluid flows through microtubes is defined as X = is/R 
where R is the radius of the tube (it is called the Knudsen number for gases, and 
in that case tg is the mean free path). Both gas and liquid flow measurements show 
dramatic enhancements in flux rates over those predicted with classical continuum 
no-slip flow models [Neto et al. 2005; Whitby and Quirke, 2007]. Majumder et al. 
[2005] observed that the measured flow rates through carbon nanotubes were four to 
five orders of magnitude greater than what was predicted by hydrodynamics mod- 
els. The authors also showed that the slip lengths, which are much greater than the 
tube diameter, are consistent with a nearly frictionless interface due to the lower 
molecular corrugation. Similar observations were reported by Cheng and Gordanio 
[2002]. 

Experiments overwhelmingly show that the presence of a slip velocity on the 
boundary can play an important role in small scale fluid devices. Theoretical solu- 
tions to various flow problems for Newtonian fluids with slip boundary conditions 
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have been obtained and are available in the literature [Karniadakis et al., 2005]. 
Fluid flows through rigid microtubes have been studied in the sHp flow regime [Sbra- 
gaglia and Prosperetti, 2007] , and this study determined that an effective Navier slip 
boundary condition is valid. Other studies have investigated steady-state [Matthews 
and Hill, 2007] and transient [Wu et at, 2008] flows through rigid microtubes with 
slip boundary conditions. Wu et at [2008] showed that the influence of boundary- 
slip on the flow behavior is different for different types of pressure gradients driving 
the flow. Since there is no wave motion in the rigid tubes, the fluids oscillate in the 
bulk, and hence the effect of wave propagation was not modeled in those studies. 

Unsteady flow characteristics and wave propagation through elastic tubes has 
been studied on the macroscale [Wormersly, 1955; Zamir, 2000]. There is no cor- 
responding study in the literature, however, for the micro- or nanoscale. Pulsatile 
microfluidics is a widely unexplored area and the difference in the dynamics at the 
macro- and microscales remains to be quantified. Pulsatile microfluidics through 
the elastic conduits is relevant to many microsystems. For example, the boundary 
tuning mechanism in surface acoustic wave driven flows (for precise control of flu- 
ids through a miniaturized network) [Chu, 2004]. Another relevant application is 
valveless pumping mechanisms. Thus, a study of fluid flows in deformable micro- 
or nanotubes should be performed to better understand some natural process such 
as respiratory mechanism in insects [Westneat et ai, 2003], flow characteristics in 
micro- and nanofluidic systems, and the efficient design of these systems. We under- 
take such a study here. In particular, we derive an analytical solution for transient 
flows through elastic microtubes with Navier slip boundary conditions. The nonlin- 
ear influence of the slip length on the material deformation of the fluid significantly 
affects the velocity profiles, fiow rates, and shear stress distributions in the system. 
Our results may be useful to researchers in the microfluidic device and biology com- 
munities for understanding the effects of pulsatile flow dynamics, elastic conduits, 
and slip velocities in internal flows at small scales. 



2. Mathematical formulation 

We consider a circular microtube with a thin elastic wall filled with a incompressible 
Newtonian fluid as illustrated in Fig.[l] The governing equations for the conservation 
of momentum of incompressible flows can be written in vector form as 




P 



1 



P 



(1) 



where fj, and p are the viscosity and density of the fluid, assumed constant for 
incompressible Newtonian fluids, u and p are the velocity and pressure fields of the 
fluid, and / represents external body forces. 
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Fig. 1. Schematic of elastic microtube with cylindrical coordinates 

These equations can be simplified according to the conditions for which we would 
like to solve them. The assumptions we make in this study are: (i) the wave speed 
of propagation is assumed to be much greater than the average flow velocity (which 
simplifies the nonlinear terms in the Navier-Stokes equation), (ii) the length of tube 
is much greater than the tube radius (which simplifies some viscous terms having 
an axial derivative), (iii) there is no external force such as gravity, (iv) the flow is 
ax:isymmetric and there is no swirling flow (no dependence on angular direction), 
(v) the thickness of the tube wall is much smaller than the tube radius (membrane 
equations for thin- walled elastic structures may be used), (vi) linearized coupling 
boundary conditions at the fluid-structure interface (which eliminates the necessity 
of a domain transformation), and (vii) the axial fluid velocity relative to the solid 
surface is proportional to the shear stress on the interface (Navier slip boundary 
conditions). These assumptions are valid in many biological and industrial processes. 

2.1. Governing equations for the fluid 

Under assumptions (i-iv) the axial and radial momentum equations, and the con- 
tinuity equation for an incompressible Newtonian fluid in cylindrical coordinates 
become [Womersley, 1955; Zamir, 2000] 

du 1 dp u. ,d^u Idu^ 
ot p oz p or^ r or 

dv 1 dp ^ fi ^d^v ^ 1 dv V ^ 
dt p dr p dr"^ r dr 



= (4) 



du dv V 
dz dr r 

where p is the fluid pressure (since only the gradient of the pressure occurs in the 
equation, it is usually called gage or transmural pressure) and u, v are the axial and 
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radial velocity components. The stress field in the incompressible Newtonian fluid 
is related to the velocity field by the following constitutive relation 

a = -pl + 2^d (5) 

where the shear stress is defined as r = 2/xd while the rate of deformation tensor d 
is related to the velocity vector by 

d^^iVu + iVuf) (6) 

where u = vcr + ucz and therefore the velocity gradient tensor becomes: 



The shear stress can then be written as 



" dv 

dr 





du " 

dr 





V 

r 





dv 





du 


L dz 


dz J 



(7) 



du dv du . 
2.2. Governing equations for the tube wall 

Using assumption (v), for thin- walled elastic circular shells the elasticity equations 
for radial and axial displacements of the tube wall are given by [Womersley, 1955; 
Zamir, 2000] 

dt^ {l-i^^)p^^dz^ Rdz' hp^ ^' 

d^V ^ Pw - Po E ri d^ 

dt^ hp^ {l-u'^)RpJ R dz' ^ ' 

where po and Pw are the values of the pressure on the outer surface and inner surface 
of the tube, respectively. Here, R and h are the reference radius and the thickness 
of the tube as shown in Fig. [l] and p^, E, and v are the density. Young's modulus, 
and Poisson ratio of the tube material. 



2.3. Fluid- structure coupling at the interface 

The dynamic forces from the fluid acting on the structure are modeled in terms of 
the shear stress in the axial direction and the pressure in the radial direction and 
give rise to the dynamic coupling conditions 

Pw-Pa=p- (11) 

where p is the transmural pressure. The kinematic boundary conditions are 

— {z,t)^u^{z,t) (13) 



December 4, 2012 2:37 WSPC/INSTRUCTION FILE ws-ijam 



6 Omer San and Anne E. Staples 

^{z,t)=v^{z,t), (14) 

where and are the axial and radial velocities of tube. Here, we linearize 
the actual boundary conditions and apply them only approximately at the neutral 
position of the wall, r = R, using assumption (vi). 

2.4. Navier slip boundary conditions 

Recent experiments and molecular simulations in microfluidics confirm that slip of 
fluid on solid surfaces occurs at small scales [Whitby and Quirke, 2007]. In this study, 
the generalized nonlinear Navier boundary condition advocated by Thompson and 
Troian [1997] is applied to the conventional continuum mechanical description of 
fluid flow considered here. The standard no-slip boundary condition is replaced by 
the nonlinear Navier slip boundary condition, wherein the slip velocity is assumed 
to be proportional to the tangential viscous stress and the degree of slip is measured 
by a non-constant slip length, tg. I.e., the axial fluid velocity relative to the tube 
wall is linearly proportional (by assumption (iv)) to the shear stress at the wall. 
The kinematic matching condition for the axial velocity can then be written using 
the Navier slip boundary condition as 

uiR,z,t) -u^{z,t) =^ -e /"-'^^^''*^ (15) 

giving 

u{R,z,t) - Uw{z,t) = -£s-Q^{R,z,t) (16) 
and the kinematic matching conditions for the radial velocity component are 

v{R,z,t)-vUz,t)=0. (17) 



3. Exact solutions 

Since the governing equations are linear we seek a solution for oscillatory flow having 
the angular frequency oj of the form [Zamir, 2000] 

p(r,.j,t) = P(r)e'"(*-^/=) (18) 
u{r,z,t) = Uir)e'^^'-'/^^ (19) 
w(r,0,t) = y(r)e"^(*-^/'=) (20) 

where c, the wave propagation speed, is a complex variable. This assumption is 
limited to cases of harmonic oscillatory motion. It can, however, be generalized in 
terms of a Fourier series expansion of diff'erent harmonics since our equations are 
linear. If we deflne the following transformation 

C = A^, A=(^)Wo, Wo=^R (21) 
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and substitute it into the governing equations we obtain the following ordinary 
differential equation for continuity 

V' + - irU = 22 

C cA 

and the following Bessel equation for axial momentum 

C^U" + cu' + c^[/ = ^ (23) 
pc 

which has a general homogenous solution of the form 

C/'^(C) = aiJo(C) + a2l^(C) (24) 

where ai, and a2 are integration constants and Jq, and Iq are the zeroth order Bessel 
functions of first and second kind [Yousif and Melka, 1997] . Since U is bounded at 
C = 0, we see that a2 = 0. The radial momentum equation becomes 

^2v" + ^V' + ie-l)V='^^^ (25) 

pRuj 

which has a homogenous solution of the form 

V'iO^b^MO + b^Y.iO (26) 

where bi, and 62 are integration constants and Ji, and Yi are the first order Bessel 
functions of the first and second kind. Similarly, V is bounded at C = 0, and 62 = 0. 
General solutions satisfying the governing differential equations then become 

P[r)=BM — ) (27) 

U{r) = ^Jo(A^) + — ^ Jo( — ) (28) 

R c[ip — IMjJ/C) c 

^ ^ , . , iujR ^ , , r . iB ^ , itor . 

V{r) = A—MA-) + — ^Ji( — ). 29 

cA H c[ip — p,u!/c) c 

If we use the following simplifications: Jo(^) ~ 1, Ji(^) ~ and p > fiui/c, 
which are valid under assumption (i), the solution of the system of differential 
equations becomes 

P(r) = B (30) 

U{r) = AMA'-) + ^ (31) 
H pc 

, ^ , , , iuiR ^ , . r , „ 

yir)=A—MA-) + B—^r. (32) 

Here, A and B, the integration constants, and c, the wave propagation speed, remain 
to be determined. It is reasonable to assume that the radial and axial oscillatory 
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movements of the tube have the same frequency as that prevaihng in the flow field, 
so that 

= Ce'"(*-"/") (33) 

77(z,i) = ZJe''"^*-^/'^). (34) 

This does not, however, imply that the wall motion in phase with oscillatory motion 
of the fluid [Zamir, 2000]. Therefore, we have additional two constant C and D that 
need to be determined. If we use the the shell equations along with the matching 
conditions given in the previous section we get a set of four homogenous hnear 
equations for A, B, C, and D with rank three which can be written as 

jjAJi(A) p. 2/1 Eh \ -iu)vE 

r, 1^ iuvE —E 

Jo(A)-|AJi(A) i ^ n } ~ ^- (^^i 

This matrix equation has a nontrivial solution if the determinant of the coefficient 
matrix is zero, and gives an equation for the wave propagation speed. After some 
algebra the quadratic dispersion relationship becomes 

ae^ +P9 + "/^0 (36) 




where 



where 



a = {l-iy^){l-g-x) (37) 
l3^g{2iy-^)-K{l-g~x)~2{l-x) (38) 
7 = g + 2«:(l-x) (39) 



Eh p^h ,4Ji(A) IJi(A) 



(1-1/2)^^^2' pji^^ i?Jo(A)'-^ AJo(A) 
Using the positive solution of Eq|36] 

^ -/? + v//32 - 4a7 



2a 



the wave propagation speed becomes 



c 



CO V (1 - ^^)^ 



(41) 



(42) 



where we use the Moens-Korteweg inviscid wave propagation speed, cq, for normal- 
ization, which is defined as 

Eh , , 
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Here, the parameter 0, and the wave propagation speed, c, are functions of the fluid 
properties and tube structural properties, geometry and slip length. We define the 
slip ratio as A = ^ and its influence on the flow characteristics is presented in this 
study. The variation of the real and imaginary parts of the wave propagation speed 
c, normalized in terms of the inviscid wave speed, cq, with frequency parameter. Wo, 
for various slip ratios is shown in Fig |ll[ Here, A = corresponds to the classical 
Womersley solution for elastic tubes. As A increases, both the real and imaginary 
parts of the wave propagation speed increase for certain range of Wo numbers. As 
the frequency increases, the imaginary part of c vanishes while the real part becomes 
the same as inviscid wave propagation speed. The effects of the wave propagation 




Wo 

Fig. 2. Variation of the real (solid lines) and imaginary (dashed lines) parts of the wave speed 
c with frequency parameter, Wo. The wave speed is normalized by the inviscid wave speed cq for 
various normalized slip ratios A = % . 



speed can be seen in the solution via the term e'"*^* If we define 

+ (44) 

C Ci C2 

where Ci and C2 are dispersion and attenuation coefficients, then the term becomes 
^ujz/c2^iui{t-z/ci) ^ Variations of the dispersion and attenuation coefficients with fre- 
quency parameter are shown in Fig[3] for different slip ratios. As the frequency 
parameter increases C2/C0 — > —00 and Ci/cq — ?> 1, and both attenuation and dis- 
persion effects vanish earlier with increasing slip ratios A. In addition to the wave 
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Fig. 3. Variation of the dispersion(ci) and attenuation(c2) coefficients witli the frequency pa- 
rameter Wo for various normalized slip ratios. Solid and dashed lines show ci/cq, and 02/ eg, 
respectively. As the frequency parameter increases both attenuation and dispersion effects vanish 
with C2/C0 — > —00 and ci/cq — > 1. 



propagation speed, as a consequence of the rank three linear system of the four 
equations, the solutions for the coefficients A, C, and D are expressed in terms of 
B as 

A-( ^ )( 2 + g(2.-l) . . . 

^-VcJo(A)^^%-2.(l-x))^ ^^'^ 

»i? g(l-g-x)-2(l-x) , 

pc^ 6'(5 - 2zy(l - x)) 
Considering the solution for pressure 

p(r,z,<) =Be*"(*-^/^) (48) 
it can be seen that the pressure gradient becomes 



dp iBuj ^,^it-zlc) 



dz c 

Prescribing the amplitude of the pressure gradient completes the analytical solutions 
derived here and closes the system of equations. To provide a reference solution, we 
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choose the ampHtude of pressure gradient to be the steady-state PoiseuiUe flow 
pressure gradient, $, setting B = We normahze our solutions with respect to 
the steady-state PoiseuiUe flow solution, which is given as 

.. = -f'5>[l-(^n (50) 
with the corresponding flow rate and shear stress 

9s = -^$, (t^). = (51) 
and we use the following centerline (r — 0) velocity to normalize the equations 

u. = -^<f- (52) 

Interestingly, the solution of the steady flow with slip boundary condition results in 
following field 

ufP^u,[l-iL)^+^-^] (53) 

9f'' = -^'i>[l + ^] (54) 

(r.)f ^ - {r^)s (55) 
where Navier slip boundary condition have been implemented. The term ^ in 



Eq. [53] represent an additional velocity associated with the Navier slip boundary 
conditions. For finite £s, as the radius of tube decreases the importance of this term 
increases. Consequently, as discussed by Wu et al. [2008], the implementation of 
slip boundary conditions adds a rigid body translation term to the axial velocity. 
The slip solution becomes a superposition of the no-slip solution and rigid body 
translation. Since rigid body translation does not lead to any change in the material 
deformation the shear stress is not affected. When slip boundary conditions are 
introduced to transient oscillatory fiow, the effect on the solution is not simply 
the addition of a rigid body translation term, and nonlinear behavior results. Wu 
et al. [2008] examined this subject by considering oscillatory flows through rigid 
tubes. In this study, we extended his study to flow through the elastic tubes, and 
by considering the axial and radial momentum equations and the cylindrical shell 
equations. If we substitute Eqs. [30p2| and the corresponding coefficients given in 
Eqs. [45]|47| into the Eq. [TSj the following normalized solutions are obtained 

'^iiA = -±[1 - G^^y-^'-/^^ (56) 

Us A2' Jo(A) J 
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I^^I^ = [G.g]e^"(*-^/^) (58) 

y^w }s 

"-^^ = - Gg'-^^V-^'-""' (59) 

Us iK^c^R ^ Ji(A) ^ ^ ' 

i? i? 

where G, the elasticity factor, an indicator of the difference between elastic and 
rigid tubes, is given as 

pcJo(A) _ 2 + g(2z.-l) 

- 0(2.(1 -X)- 5) ^''^ 
where comes from dispersion relationship, Eq. |36[ and is a fmiction of the slip 
parameter x- It can clearly be seen that the effect of slip boundary conditions is 
nonlinear in this case. The variation of the real and imaginary parts of the elasticity 
factor with Wo are illustrated in Fig.|4]for different A values. As a special case, when 
X = 0, corresponding to no-slip boundary conditions, the solution reduces to the 
classical Womersley's solution for elastic tubes. Additionally, for G = 1 and c — oo 
the solution reduces to Womersley's solution for the rigid case. 



4. Results 

Solutions of the governing equations for oscillatory flow through the elastic tubes 
were derived using Navier slip boundary conditions in Section 3, and the correspond- 
ing analytical solutions given in Eqs. [56]|6T In the following results, the constant 



values Co = lOm/s and k = 0.1 are used. The parametric variables are Womersley 
number. Wo, and the slip ratio, A. 

The influence of the slip length on the slip velocity at the surface is shown 
in Fig. [5] for various Womersley numbers. Results for the rigid wall case are also 
included. In order to make the presentation clear we present four different figures 
with different A scales. As can be seen in the figures, the slip velocities at the wall 
increase rapidly with increasing slip length and then, after a certain critical value 
of the slip length, reduce slowly. Comparing the elastic and rigid tube results, one 
can see that the maximum slip velocity can be obtained at a smaller slip ratio for 
the elastic tubes. It can also easily be seen that for large Wo there is a contribution 
to the slip velocity, which is due to axial deformation. For small Wo we extend the 
range of A to be able to see the full behavior of the slip velocity. This can be related 
to nanoscale experimental observations that the slip length can be much bigger than 
the tube radius at the nanoscale. Our results demonstrate the same conclusions. We 
note, however, that these large values of A may have physical meaning only for liquid 
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Wo 

Fig. 4. Variation of the real (solid lines) and imaginary (dashed lines) parts of the elasticity factor 
G for different values of A. 

flows at the nanoscale, which correspond to nearly frictionless flows. Majumder et 
al. [2005] measured the shp length of a carbon nanotube with a diameter of 7nm 
to be about 3 — 70/xm for different liquid flows (water, ethanol, propanol, hexane, 
and decane). In the rest of the paper we restrict the results to those for small A to 
make the analysis applicable for a wide range of systems. 

The amplitude of the percent difference between the peak flow rate in an elastic 
tube and that in a rigid tube is shown in Fig. [6] for various slip ratios. By increasing 
the slip length at the wall the difference between the rigid and elastic tubes reaches 
its peak values at a smaller value of Wo. We can also interpret this as being for a 
fixed frequency and smaller size microtubes (Wo — \/ReSt, where Re and St are the 
Reynolds and Strouhal numbers) . It is already established that the wall movements 
in an elastic tube make it easier for the flow to move through the tube than through 
a rigid tube [Zamir, 2000]. Here we show, in addition, that the the velocity slip at 
the surface also has the same behavior, with increasing sensitivity to Wo. Therefore, 
the driving frequency has even more of an effect in controlling the flows through 
smaller devices in pulsatile microfluidics. 

The normalized shear stress amplitudes at the wall are shown in Fig. [7] for both 
elastic and rigid walls. The solid and dashed lines show the slip (A = 0.1) and no-slip 
(A — 0) cases, respectively. It can be seen that the shear stress at the walls in elastic 
tubes is greater than that in rigid tubes for smaller Wo. There is always a critical 
Wo at which the rigid and elastic solutions intersect, and this critical Wo number 
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10 r 250 r 




Fig. 5. Influence of the slip ratio A on the normahzed sHp velocity at the wall. Solid lines show 
the elastic wall case and dashed lines show the rigid case: (a) Wo = 10, Wo = 7.5, Wo = 5; (b) 
Wo = 4, Wo = 3, Wo = 2; (c) Wo = 1, Wo = 0.75, Wo = 0.5; (d) Wo = 0.3, Wo = 0.2, Wo = 0.1. 



decreases with increasing slip ratio. It can be also seen that the shear stress at the 
wall decreases with decreasing Wo and nearly coincides with the frictionless flow 
measurements at the smallest values. 

In the previous figures we analyzed the amplitude of the flow rate, the slip- 
velocity, and the shear stress at the walls, and compared the results for both elastic 
and rigid tubes with different slip ratios. For those cases, the exponential terms in 
the solutions took unit value (i.e., z = 0, t = 0). Now, we introduce the axial tube 
length, z : [0, Lj], and the time with respect to one oscillatory period, t : [0,T]. 
We assume that the the tube length is much greater than the tube radius, and 
that the wavelength, L^, is also much greater than the tube radius {Lt >> R). 
Therefore, in the limit of the long-wavelength assumption, the wavelength i^, can 
be of the same order or greater than the tube length Lt . The wavelength is defined 
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as Lm = cT = 2ttc/uj. 

The particle trajectories of the tube waU using radial and axial displacements 
given in Eqs. 60 and 61 are shown in Fig. [8] for Wo = 1 and Wo = 3. The dis- 
placements are normalized by (i-e., the axes in the figure are ^/[ii?<i>/(pwc)], 
ri/[iR^/{pujc)]). As shown in the figure, the trajectories are elliptic and their areas 
increase for increasing slip ratios A. The solid line shows the no-slip (A = 0) case, 
and the dashed and dash-dotted lines show the slip cases for A = 0.1 and A = 0.2. It 
can be seen that the amplitudes of the displacements increase with increasing slip 
ratio. Additionally, the aspect ratios of the elliptical trajectories become smaller 
with increasing Womersley numbers. The orientation of the elliptical axes is also 
changing by increasing the slip ratio. For higher slip ratios (i.e., smaller scales), we 
can see that the slope of the longitudinal axis of the elliptical trajectories gradually 
decreases. 

The following results for the axial velocity are presented to characterize the fiow 
properties in elastic microtubes. The normalized axial velocity profiles according to 
Eq. [56]for different frequency parameters and tube lengths are shown in Figs. [9p2| 
We include the solutions corresponding to the rigid case as well. We also compare 
with the no-slip and Navier slip boundary conditions cases. In these figures, the 
slip ratio is chosen to be A = 0.1, so that the slip length at the wall is ten times 
smaller than the tube radius. In each subfigure the line styles and colors indicate 
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Fig. 7. Normalized wall shear stresses amplitudes. The solid and dashed lines show the elastic 
and rigid wall cases, respectively. 
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Fig. 8. Parti cle traj ecto ries of the tube wall using both radial and axial displacements. These are 
given by Eqs. 60 and 61 which are normalized here by J*-, 



the following cases: the red dash-dotted line: no-slip elastic wall, the black dash- 
doubledotted line: no-slip rigid wall, the blue solid line: slip elastic wall, and the 
green dashed line: slip rigid wall. The profiles are presented at two discrete axial 
locations and for two different phases. 

In Fig. [9] and Fig. [TOj Wo = 1 and Wo = 3, respectively, for the tube length 
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Fig. 9. Normalized velocity profiles for Wo = 1 and Lt/L^ = 1 at (a) t/T = 0, z/Lt = 0; (b) 
t/T = 0, z/Lt = 1/2; (c) t/T = 1/4, z/Lt = 0; and (d) t/T = 1/4, z/Lt = 1/2. Red dash-dotted 
line: no-slip elastic wall, black dash-double dotted line: no-slip rigid wall, blue solid line: slip elastic 
wall, green dashed line: slip rigid wall. 




Fig. 10. Normalized velocity profiles for Wo = 3 and Lt/L^i = 1 at (a) t/T = 0, z/Lt = 0; (b) 
t/T = 0, z/Lt = 1/2; (c) t/T = 1/4, z/Lt = 0; and (d) t/T = 1/4, z/Lt = 1/2. Red dash-dotted 
line: no-slip elastic wall, black dash-double dotted line: no-slip rigid wall, blue solid line: slip elastic 
wall, green dashed line: slip rigid wall. 
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Fig. 11. Normalized velocity profiles for Wo = 1 and Lt/L^u = 0.1 at (a) t/T = 0, z/Lt = 0; (b) 
t/T = 0, z/Lt = 1/2; (c) t/T = 1/4, z/Lt = 0; and (d) t/T = 1/4, z/Lt = 1/2. Red dash-dotted 
line: no-slip elastic wall, black dash-double dotted line: no-slip rigid wall, blue solid line: slip elastic 
wall, green dashed line: slip rigid wall. 




Fig. 12. Normalized velocity profiles for Wo = 3 and Lt/L^ = 0.1 at (a) t/T = 0, z/Lt = 0; (b) 
t/T = 0, z/Lt = 1/2; (c) t/T = 1/4, z/Lt = 0; and (d) t/T = 1/4, z/Lt = 1/2. Red dash-dotted 
line: no-slip elastic wall, black dash-double dotted line: no-slip rigid wall, blue solid line: slip elastic 
wall, green dashed line: slip rigid wall. 
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to wavelength ratio = 1. When Wo is lower as in Fig. joj the velocity profiles 
become close to parabolic profiles. As discussed earlier, the presence of slip at the 
solid surface has only translational effects for constant pressure driven flows. When 
we consider pulsatile flows, however, the slip length has a nonlinear effect. Nonlinear 
effects of the slip length are obvious for the elastic wall case as well, as is demon- 
strated in the figures. We see that the solutions for a rigid tube do not have any 
wave motion, so the velocity profiles at different sections are always the same at the 
same time (the wave propagation speed, c — oo). We can see the wave motion in 
the elastic tubes, however, due to the finite wave propagation speed; this phenom- 
ena can easily be seen because of the comparable tube lengths and wavelengths. 
The corresponding velocity profiles for the same tube length ratio but moderate 
frequency parameter Wo = 3 are shown in Fig. |10| Since the frequency parameter 
is moderate we see that the velocity profiles are not parabolic and the amplitude is 
reduced at the centerline. Furthermore, if we compare the slip and no-slip solutions, 
we see that this reduced centerline velocity effect increases with slip velocity at the 
wall. This is due to fluid inertial effects. For larger Wo the flow has inertia acting in 
the opposite direction when the pressure gradient is reversed. Therefore, it will take 
some time before the pressure gradient can reverse the direction of the inertia. This 
time delay introduces a phase shift between the fluid motion and the pressure gra- 
dient. In a confined domain such as a tube, smaller velocities with correspondingly 
lower inertia towards the wall have smaller phase shifts than the fiuid velocities in 
the center of the tube. 

Fig. [TT] and Fig. [T2| show the velocity profile for Wo = 1 and Wo = 3 for the tube 
length ratio of =0.1. Since the wave lengths are greater than the tube length, the 
waves complete their cycles more quickly in space, and the phase difference between 
the rigid and elastic solutions decreases. For the rigid tub case, the waves complete 
their cycle at an almost infinite speed, and so the phase difference between the rigid 
and elastic tube solutions becomes less important. It can also be seen that for small 
Wo the fluid reaches a steady state at all times before the pressure changes, and 
therefore a quasi-steady state is reached, and, correspondingly, a Poiseuille velocity 
proflle. 

Finally, we consider the normalized axial velocity contours for elastic microtubes 



with two different tube length to wave length ratios. Fig. 13 and Fig. 14 present the 
Lt/ Lw = 1 and Lt/ = 0.1 cases, respectively. These figures include axial velocity 
fields for two different frequency parameters. Wo = 1, and Wo = 3, and a slip ratio 
of A = 0.1 is used in each of them. Concerning the harmonic oscillations, we plot 
only one phase state of the axial velocity fields of interest in one time period (t = 0) . 



As seen in Fig. [T3j the wave completes its cycle in space if the wavelength and tube 
length are equal. As discussed earlier, the velocity profiles become more parabolic 
when Wo is small. Increasing Wo affects the flow fleld by effectively increasing the 
ratio of oscillatory inertial forces to shear forces. The effect of the slip length as a 
function of Wo can be seen in the figures, and it is clear that it interacts nonlinearly 
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Fig. 13. Normalized axial velocity contours inside a tube with Lt/L^ = 1 for (a) Wo = 1, and 
(b) Wo = 3. 




Fig. 14. Normalized axial velocity contours inside a tube with LtjL^ = 0.1 for (a) Wo = 1, and 
(b) Wo = 3. 
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with the pulsatile flow dynamics. 



5. Conclusions 

The influence of a slip velocity on the flow characteristics in circular clastic mi- 
crotubes is investigated using a continuum approach with Navier slip boundary 
conditions. Analytical solutions of the pulsatile flow dynamics in a small tube are 
found by taking into account the coupling between the fluid motion and the elastic 
deformation of the tube wall. The formulation derived here can be reduced to the 
Womersley solution for elastic tubes as a special case by setting the slip length 
equal to zero. In addition to that, by setting the elasticity factor G — 1 along 
with a infinite wave propagation speed, we also recover the classical solution for 
oscillatory flows through rigid pipes. The influence of slip boundary conditions on 
the flow through tubes for pulsatile flows is qualitatively different than for steady 
flow with a constant pressure gradient, wherein slip boundary conditions only result 
in an added translational effect, without changing the material deformation of the 
fluid. With an unsteady pressure gradient, however, the influence of slip boundary 
conditions becomes nonlinear and affects the flow rate, velocity proflle, and shear 
stress, by altering the material deformation in the fluid. 

There are large differences between flows in elastic and rigid tubes, especially 
if the tube length and the wavelength are on the same order. It is already well 
established that the wall movements in an elastic tube make it easier for the; flow 
to move through the tube than through a comparable rigid tube. Solutions for rigid 
tubes do not have any wave motion, so velocity profiles at different sections are 
the same for the same point in the oscillation period (the wave propagation speed 
c — > oo). There arc wave motions in clastic tubes due to finite wave propagation 
speeds, however, and this phenomena can be easily seen when the tube length is 
comparable to the wave length. 

In this study, we demonstrated that the presence of slip at the elastic microtube 
wall enhances flow through the tube. This effect is not just an additional transla- 
tional velocity, but shows nonlinear behavior and affects the material deformation 
in the flow fleld. Additionally, we confirmed that the elasticity of the tube wall 
also has a positive effect on the flow rate. There is always an increase in flow rate 
for elastic tubes compared to the inelastic case, and the peak value of the percent 
change in flow rate increases with increasing slip length. The effective range of the 
frequency parameter Wo decreases, however, with increasing slip length, and the 
corresponding value of Wo at this peak value also decreases with increasing slip 
length. In summary, the presence of velocity slip at the microtube surface has im- 
portant, nontrivial effects on the velocity profiles. It changes the amount of the 
phase shift between the fluid inertial motion and the pressure gradient, and alters 
the material deformation of the media. 
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